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ABSTRACT 



Aims. We study the diffuse X-ray emission observed in the field of view of the pulsar B 0540-69 in the Large Magellanic Cloud 
(LMC) by XMM-Newton. We want to understand the nature of this soft diffuse emission, which coincides with the superbubble in 
the Hn region N 158, and improve our understanding of the evolution of superbubbles. 

Methods. We analyse the XMM-Newton spectra of the diffuse emission. Using the parameters obtained from the spectral fit, we 
perform calculations of the evolution of the superbubble. The mass loss and energy input rates are based on the initial mass function 
(IMF) of the observed OB association inside the superbubble. 

Results. The analysis of the spectra shows that the soft X-ray emission arises from hot shocked gas surrounded by a thin shell of 
cooler, ionised gas. We show that the stellar winds alone cannot account for the energy inside the superbubble, but the energy release 
of 2 - 3 supernova explosions in the past ~1 Myr provides a possible explanation. 

Conclusions. The combination of high sensitivity X-ray data, allowing spectral analysis, and analytical models for superbubbles bears 
the potential to reveal the evolutionary state of interstellar bubbles, if the stellar content is known. 
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1. Introduction 



Early observations in the radio and the optical have shown that 
the interstellar medium (ISM) in the Milky Way mainly consists 
of cool clouds (T ^ 10 2 K) of neutral hydrogen embedded in 
warm (T 10 4 K) intercloud medium of partially ionised hy- 
drogen. Since the 1970s, observations in the ultraviolet (UV) and 
X-rays showed the presence of hot gas at coronal temperatures 
(T a 10 5-6 K) in the ISM. The heat source of the ISM are mas- 
sive OB stars, which inject energy through their radiation, stellar 
winds, and finally by supernova explosions. As these processes 
are often correlated in space and time, superbubbles with sizes 
of typically 100 - 1000 pc are created in the ISM. Therefore, 
supernova remnants (SNRs) and superbubbles are among the 
prime sources that control the morphology and the evolution of 
the ISM, and their observation is of key interest to understand 
the galactic matter cycle. However, they radiate copiously in the 
soft X-rays below 2 keV, an energy range that is difficult to study 
in the Milky Way because of absorption by the Galactic disk. 

The LMC, which is a dwarf irregular, but shows indica- 
tions for spiral structures, is one of the closest ne ighbours of 
our G alaxy. Its proximity with a distance of 48 kpc (iMacri et alJ 
2006) and modest extinction in the line of sight (average Galactic 
foreground 7V H = 1.6 x 10 21 crrT 2 ) make it the ideal laboratory 
for exploring the global structure of the ISM in a galaxy. The 
well-known and best studied extended emission region in the 
LMC is the 30 Doradus region and the region south of it, which 
harbor star formation sites, superbubbles, and SNRs. ROSAT 
data of the superbubb l es in the LMC have been stud ied in de- 
tail bv.e.g JChu etal] (119951) and lDunne etai1(l200ll) . 



N 158 (Henize 1956) is one of the H n regions in these active 
regions of the LMC. It is elongated in the north-south direction 
and consists of a superbubble in the north and a more compact 
bright region in the south . It is known to host tw o OB associa- 
tions LH 101 and LH 104 dLucke & Hodgel 19701) . While LH 101 
in the southern part of N 158 seems to power the very bright re- 
gion in Ha, LH 104 is found in the superbubb le in the north- 
ern p art of N 158 and is dominated by B stars dSchild & Teston 
Il992l) . main ly consisting of a young population w ith an age of 
2-6 Myr dTestor & Niemelalll998h . iDunne et all (1200 ll) have 
analyzed the ROSAT data and suggested that the X-ray emis- 
sion seen at the position of N 158 is associated with the Hn 
region. N 158 is located near the X-ray bright pulsar (PSR) in 
the LMC B 0540-69, which has been observed for calibration 
purposes for the X-ray Multi-Mirror Miss ion XMM-Newton 
ilJansen et alJl200U lAschenbach et all 12000). The field of view 
of the European Photon I maging Cameras (EPICs. IStrtider et"a"fl 
200 1 : tTurner et al.1 1200 1 ) of these observations when performed 
in Full Frame mode, covers the northern part of N 1 58 and allows 
to study the X-ray emission from the superbubble. 



2. Data 

The pulsar B 0540-69 in the LMC is a Crab-like pulsar 
with a pulsar wind nebula (PWN), which has been spatially 
resolved and st udied with the Chandra X-ray Observatory 
dPetre et al.ll2007l) . To study the diffuse emission in the vicinity 
of B 0540-69, we have chosen those observations for which the 
EPICs were operated in Full Frame mode. The observation IDs 
are 0117510201, 0117730501, and 0125120101. The observa- 
tions were all carried out using the medium filter. Starting from 
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Table 1. XMM-Newton data used for the analysis. All the anal- 
ysed data were obtained in Full Frame mode using the medium 
filter. 



Obs. ID 


EPIC 


Start Date 


Effective Exposure 
[ksec] 


01175102 


PN 


2000-02-11 


8.3 


01175102 


MOSl,2 




3.5 


01177305 


PN 


2000-02-17 


8.3 


01177305 


MOSl,2 




9.8 


01251201 


PN 


2000-05-26 


29. 


01251201 


MOSl,2 




27. 



the observational data files (ODFs), the data are processed with 
the XMM-Newton Science Analysis System (XMMSAS) ver- 
sion 10.0.0. For EPIC PN, only single and double pattern events 
are used, whereas for the MOS 1 and 2, singles to quadruples are 
selected. The exposure times that we obtain after filtering out the 
time intervals with background flares are listed in TableQ] 



2.1. EPIC image 

After filtering out the background flares, we created a mosaic 
image out of the Full Frame mode data of EPIC PN, MOS 1 , and 
MOS2 for all three observations (Fig.[TJ. To enhance the not so 
bright diffuse emission, we filtered out all point sources that have 
been found in a source detection routine as well as the so called 
out-of-time events of EPIC PN. The images have been smoothed 
using a Gaussian filter. The mosaic image shown in Figure [TJ is 
a true color image using the colors red for the 0.3 - 0.8 keV 
band, green for 0.8 - 1.5 keV, and blue for 1.5 - 2.3 keV. The 
extended emission of the interstellar gas is clearly soft with no 
emission above ~3 keV. The comparison with the Ha image in 
Figure [2] shows that the relatively bright extended region in the 
south of the P SR coincides well with a superbubble in the H n 
region N 158 ( Henizdl 19561). w hich contains the OB association 
LH 104 dLucke & Hodgelll97dl) . In order to study the spectral 
properties of the diffuse emission, we have selected two regions: 
region 1, which covers the brighter spot in the east of the PSR, 
and a region that covers the superbubble in the H n region N 158. 
The regions are shown in the left panel in Figure [2] The PSR 
and the PWN with an extent of about 1 ' is completely removed 
from the data. The soft, extended emission east to the PSR is not 
directly connected to the PWN and has, as we will see in Section 
12.2.41 a perfectly thermal spectrum. Therefore, we assume that it 
is not related to the PWN. 



2.2. EPIC spectra 

For the spectral anaylsis of an extended diffuse emission the con- 
tribution of the background is significant. As the emission fills 
a large part of the detector, we are not able to extract a local 
background close to the source emission. We have to consider 
the following: As the effective area of the mirrors depends on 
the off-axis angle, photons are subject to vignetting while parti- 
cles are not. The high-energy particles that interact with material 
surrounding the detector, however, produce fluorescence, which 
varies with position on the detector, especially for the PN detec- 
tor. In addition, the spectral response depends on the position on 
the detector. A deta iled description of the X M M-Newton back- 
ground is given by Read & Ponman (2003) and ICarter & Read! 



(120071) . and a comparison of the diff erent methods to e stimate 
the background can be found in, e.g. JSasaki et al.1 (120041) . 



2.2.1. Blank sky background 

The XMM-Newton EPIC Background working group has cre- 
ated the s o-called blank sky dat a for each EPIC and CCD read 
out mode (ICarter & Readll2007h . The blank sky data have been 
merged from data of different pointings after the point sources 
were eliminated from the data. These data sets comprise the de- 
tector background and an average cosmic X-ray background. 

Before extracting the spectra, we have first corrected both 
the observed data and the blank sky data for vignetting using 
the XMMSAS command evigweight. The auxiliary response 
file (ARF) and the response matrix file (RMF) are then created 
assuming that the source is on-axis. The background spectrum 
is extracted from the blank sky data at the same location on the 
detector as the source spectrum. 



2.2.2. Background from the same data 

For comparison, we also extracted a background spectrum from 
the same data as the source spectrum, using a region in the dark 
part north of B 0540-69. However, after the subtraction of the 
background, the spectrum of the soft diffuse emission still has 
a hard tail and is over-corrected exactly at the energies of the 
fluorescence lines. Therefore, as expected, extracting the back- 
ground from the same data at a different position on the detector 
seems to be inappropriate. 



2.2.3. Closed filter wheel data 

Another way to deal with the background is to use a local back- 
ground from the same obervation, but take care of the detector 
background by using the closed filter wheel data, also supplied 
by the XMM-Newton EPIC Background working group. For es- 
timating the X-ray background, i.e., additional emission that is 
typical for the observed area and might also contribute to the 
spectrum of the superbubble, we extract a region next the super- 
bubble in the east, which shows faint diffuse emission. After sub- 
tracting the closed filter wheel spectrum extracted for each re- 
gion at the corresponding position of the detector with the same 
shape, the spectrum of the faint diffuse emission in the east is fit- 
ted with a thermal model. This fitted spectrum is then included 
in the model of the spectra of the superbubble and region 1 . 

To verify if the emission in the east is suitable for the use 
as the local X-ray background, we estimate the flux for the su- 
perbubble, for the region east to it used as the local background, 
and for the blank sky background. We assume a plane-parallel 
shock model for the superbubble and for the eastern region a 
non-equilibrium ionisation (NEI) model, both with LMC abun- 
dances (see Sect. l2.2.4l for details). A combination of a thermal 
and a non-thermal spectrum is assumed for the blank sky back- 
ground. From these spectra, we get the following fluxes with 
90% confidence errors in brackets: Fo2-3 0kev(superbubble) = 
1 .4 (1 .2 - 1 .5) x 10~ 12 erg s" 1 cirT 2 , F .2-3.o kev(east) = 4.8 (4.0 - 
5.3) x 10~ 13 — — " 2 
1.6) x 10 

the blank sky data is about 10% of that of the superbubble. The 
emission from the ISM in the east is about three times higher 
than the blank sky. If we use the blank sky data for estimating 
the X-ray background, we might underestimate the background. 
However, if we use the local emission, we will overestimate the 



erg s cm 2 , Fo.2-3.0 k e v(blank sky) = 1.5 (1.4 - 
13 erg s crrT 2 . The flux of the X-ray background in 
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Fig. 1. XMM-Newton EPIC mosaic image of the PSR B 0540-69 and its surroundings in true color presentation (red: 0.3-0.8 keV, 
green: 0.8 - 1.5 keV, blue: 1.5 - 2.3 keV). The bright X-ray emission from PSR B 0540-69 at ~RA = 05 ft 40 m , Dec = -69° 20', 
other point sources, and the out-of-time events have been removed from the data. The arc-shaped features in the south are caused by 
stray light from the bright X-ray source LMC X- 1 . The position of the superbubble in the H n region N 1 58 is shown with a dashed 
line. 



background as we may expect hot gas also in the ISM next to 
the superbubble because of possible breakouts. The difference 
between the flux of the eastern region and the blank sky is about 
20% of the superbubble flux, which needs to be taken into con- 
sideration as an additional uncertainty for the flux of the super- 
bubble. 

Otherwise, the results for the spectral fits obtained with the 
closed filter wheel data and those obtained with the blank sky 
data are consistent with each other within the confidence range 
of the spectral fit parameters. Therefore, in the following, we 
will discuss the fit results obtained with the blank sky data. 

2.2.4. Spectral fits 

The spectra are fitted with thermal plasma models in XSPEC. 
We fitted the spectra of the different EPICs simultaneously with 
model parameters that are linked to each other. In both spec- 
tra emission peaks are found in the energy interval between 
0.5 keV and 1.5 keV which can be interpreted as emission 
lines of highly ionised elements. The X-ray emission is ab- 
sorbed by the Galaxy in foreground (Nu fg) and the matter in 
the LMC in the line of sight (Nhmc) w ith abundances half 
of solar values dRussell & Dopital |1992|) . The Galactic fore- 
ground column density is N ufe = 7.0 x 10 20 crrT 2 in this region 
(Dick ey & LockmarJ 1 19901). The s pectra are better reproduced 
by a NEI model dBorko wski 2000) than by the spectral mod- 
els assuming collisional ionisation equilibrium (CIE). Figure [3] 



Table 2. Spectral parameters obtained from the fits of the EPIC 
data for region 1 and the superbubble. The numbers in brackets 
are 90% confidence ranges. The parameter norm is the normali- 
sation of the model in XSPEC and corresponds to the emission 
measure. For the other parameters, see Section l2.2.4l 



Parameter 




Region 1 


Superbubble 


NEI 


Whlmc [10 22 


cm J 


0.14 (0.08- 


- 0.23) 


0.43 (0.38 - 


-0.51) 


kT [keV] 




0.54 (0.51 - 


-0.58) 


0.91 (0.74 


-1.1) 


t [10" scrrT 


3 i 


2.0(1.6- 


2.4) 


0.14(0.12- 


-0.18) 


norm (1CT 4 ) 




2.0(1.6- 


2.4) 


9.5 (8.0- 


13.) 


//d.o.f 




94.81/61 = 


1.55 


185.06/106 


= 1.75 


NPSHOCK 


Whlmc HO 22 


cm 2 ] 


0.19 (0.12- 


-0.32) 


0.29 (0.20 - 


-0.33) 


kT a [keV] 




0.55 (0.00 - 


-0.59) 


1.0 (0.74- 


-1.2) 


kT b [keV] 




0.51 (0.00- 


- 0.59) 


0.13(0.00- 


- 0.57) 


T[10 n s cm- 


3 ] 


8.8(5.4- 


12.) 


1.3 (1.0- 


2.3) 


norm (10~ 4 ) 




5.0 (4.2- 


6.3) 


11. (9.- 


13.) 


* 2 /d-o.f 




91.02/60 = 


1.52 


131.79/99 


= 1.33 



shows the MOS1 spectra of the Obs. ID 01251201 with the best 
fit models. In the following, 90% errors are given for the param- 
eters. 



4 Manami Sasaki et al.: XMM-Newton observations of N 158 




Fig. 2. A zoom in on the XMM-Newton EPIC mosaic image with regions used for the spectral analysis (black dashed line) and Ha 
contours (MCELS, left) and the Ha image with the same contours (right). 



Since the first fits with CIE yielded no satisfactory results, 
we use the NEI model in XSPEC with the effective temperature 
kT, the ionisation timescale r = n e t, and the abundance as fit pa- 
rameters. The ionisation timescale t is an indicator for the state 
of the plasma after the gas has been ionised: For small ionisa- 
tion timescales t = n e t < 10 12 s cirT 3 , the electrons and ions are 
still not in thermal equilibrium. Also for this model component, 
we use the sub-solar LMC abundance £lmc = 0.5. The fit of the 
spectra of the superbubble emission is still not very good with 
a reduced^- 2 of 1.75. The emission is most likely caused by the 
shock from the stellar winds in the superbubble. Therefore, we 
also fit the spectra using a plane-parallel shock plasma model 
with separate ion and electron temperatures, kT a and kTb, re- 
spectively rNPSHOCK. lBorkowski et al.ll200ll) . The fit parame- 
ters are given in Table [2] 

As can be seen in Table [2] kT a and kTb for region 1 are al- 
most equal and the^ 2 values for the NEI and the NPSHOCK fits 
are comparable. Therefore, the usage of the NPSHOCK model 
does not seem to be necessary for region 1 . In addition, the value 
for r is > 10 11 s cnr 3 for both NEI and NPSHOCK models, in- 
dicating that the gas in region 1 seems to be closer to CIE than in 
the superbubble in N 158. The temperature of kT = 0.54 keV de- 
termined for region 1 is higher and more accurate t han the result 
of the ROSAT PSPC analysis for this field (0.3 keV. ISasaki etaD 
|2002|) . This shows that we were only able to see the overall char- 
acteristic of a larger region with PSPC, whereas with EPIC we 
can now resolve smaller (few arcminutes) structures of the hot 
ISM. 

The hot plasma in the superbubble of the Hn region N 158 
seems not to be consistent with thermal equilibrium, indicated 
by the low ionisation timescale r = 1.4 x 10 10 s cirT 3 of the 
NEI fit. The temperature of the NEI fit and the ion temperature 
of the NPSHOCK fit are relatively high: kT = 0.91 keV and 
kT a — 1 .0 keV, respectively, while the electron temperature is 
one order of magnitude lower: kTb = 0.13 keV. Additionally 
to the Galactic foreground column density, a relatively high ab- 
sorption column density of A^hlmc.niss = 4.3 or 2.9 x 10 21 cm -2 
was determined. In comparision, the tot al column density in the 
LMC is JVhlmc = 1.0 - 5.5 x 10 21 cnr 2 (B rims et all2005l) . This 



corroborates that the diffuse X-ray emission most likely arises 
from inside the shell of the superbubble. 



3. Discussion 

While in the east of the PSRB 0540-69 the diffuse emission 
seems to arise from ionised gas close to thermal equilibrium 
with a temperature comparable to the value determined from 
the ROSAT PSPC spectrum, in the south a smaller ionisation 
timescale and higher temperatures are determined. Since the 
emission coincides spatially with the superbubble in the Hn 
region N 158, we conclude that the origin of the diffuse X- 
ray emission is the hot gas within the interstellar bubble. The 
gas in the bubble interior is shocked by stellar winds, and the 
cooler outer rim is visible as an Hn region. The total unab- 
sorbed X-ray luminosity of the bubble is Lx(0.2 - 10.0 keV) 
= 1.5 xlO 36 erg s _1 . Stars with masses over 25 M are lumi- 
nous X-ray emitters with luminosities of Lx(0.2 - 10.0 keV) 
» 10 33 erg s" 1 . In LH 104 there are 16 stars with masses above 
25 M dTestor & Niemelalll998l see also SectETHfor details), 
which account for L x (0.2 - 10.0 keV) * 1 - 2 x 10 34 erg s -1 . 
This is two magnitudes lower than the emission from the bubble 
and can be neglected in the following discussion. 

The comparison of the X-ray emission with the Ha shell in 
Figure [2] shows that there is additional X-ray emission outside 
the Hff shell which might indicate that some hot g as escaped 
the su perbubble. This has also been suggested by iDunne et al.l 
(2001) based on the analysis of the ROSAT data, which had al- 
ready shown that the X-ray emission is not confined by the Ha 
shell. However, we are not able to rule out that some projection 
effect might make the Ha look smaller than the extent of the 
X-ray emission. 

3. 1 . Results from the spectral analysis 

For further studies, we use the NEI fit results for region 1 and 
NPSHOCK fit results for the superbubble. The emitting volume 
can be approximated by an ellipsoid, although it is deformed 
in the south. We derive the radii from the EPIC mosaic image: 
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Fig. 3. Spectra of the diffuse emission in the field of view of PSR B 0540-69 extracted from the XMM-Newton EPIC MOS 1 data 
of the Obs. ID 01251201. The left panels show the spectra of region 1 located east of the PSR, the right panels show the spectra of 
the superbubble in the Hn region N 158. 



a = 4' ± 1' = (56 + 14) pc, b = 3' ± 1' = (42 + 14) pc 
(Dlmc = 48 kpc). We are not able to detemine the third ra- 
dius of the assumed ellipsoid. Therefore, we assume that the 
superbubble is oblate and has a configuration like a disk per- 
pendicular to the plane of the sky, i.e., also perpendicular to the 
disk of the LMC (c = a). The volume of the bubble is then 
V = 4/3 n abc = (1.62 ± 0.03) x 10 61 cm 3 . With the LMC metal- 
licity £lmc = 0.5, n e = (1.2 + 0.013£lmc) « ~ 1.21 «, with n 
being the hydrogen density. Therefore, the normalisation of the 
spectral fit is 



1 



10 14 x4^ MC 



10^x4^ MC 



ndf V 



1.21 n 2 fV 



= 4.4 x 10" 62 n 2 / V [cm" 5 ] (1) 
and the gas density within the bubble can be estimated as: 



n = 4.8 x 10 30 



norm r 
X „/— ttt- [cm J ]. 



(2) 



With norm = (1.1 ± 0.2) x 10~ 3 from the fit with the NPSHOCK 
model, we get n = (4.0 + 0.7) x/~ I/2 x 10~ 2 cirT 3 . If / < 1 then 
the density of n — 4.0 x 10~ 2 cirT 3 is a lower limit. However, 
the angular resolution of the X-ray data does not allow us to un- 
ambiguously determine the filling factor. For a young interstellar 
bubble like this particular case in N 158 the filling factor can be 
assumed to be / « 1 . 

With temperature T a = (1.0 + 0.2) keV and density n given, 
the pressure of the gas is: 

P/k = (n e + l.ln)f- l/2 T a = 231nf- l/2 T a 
= (1.1 +0.3)x/- 1/2 



x 10 6 cm" 3 K. 



(3) 



While the pressure of the Galactic ISM is thought to be P/k = 
10 3 ~ 4 cm -3 K, star forming regions in general have higher pres- 
sures of the order of P/k = 10 5 " 6 cm" 3 K. For a galaxy l i ke the 
LMC with high star formation rate. lOev & G arcfa-Segura (2004) 
estimated an ISM pressure of P/k * 10 5 cm" 3 K. Thus the pres- 
sure inside the superbubble in N 15 8 is about ten times higher 
than i n the surrounding hot ISM. de Avil lez & Brei tschwerdt 
(2005) have performed a 3-D simulation of the ISM including 
the effect of magnetic fields and obtained a map of the distri- 
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bution of temperature, pressure, magnetic field, etc. They have 
shown that in regions where the temperature is about 10 K, 
the pressure is P/k = 10 4 ~ 5 cm" 3 K, while it can reach P/k > 
10 5 cirT 3 K in the interior of hot bubbles. This is in agreement 
with the pressure that we obtain from the X-ray spectrum of the 
superbubble. 

The part of the LMC in which N 158 is located in general 
shows faint diffuse X-ray emission indicative of hot ISM. The re- 
gion which we call region 1 is particularly bright and allows us to 
estimate the foreground column density that causes absorption of 
the soft X-rays. The shell of cooler gas around the superbubble 
should make an additional absorbing component. Using the ab- 
sorbing column density determined for the region 1 (besides the 
Galactic column density) as the mean LMC value A^HLMC.regioni 
= (1.4 + 0.9) x 10 21 cm" 2 , the column density of the shell 
around the superbubble is A^h shell = Nhlmcniss-A'h LMC, region 1 — 
(1.5 + 1.3)xl0 21 cnr 2 . 

3.2. Analytic estimates 

With the values derived from spectral fitting for pressure and 
density inside the bubble, we infer a thermal energy content of 
E(t) = 3/2xPxfxV = 3.6x/ I/2 xl0 51 erg. We estimate a mean 
shell thickness of ~0'5, i.e., ~ 10% of a, corresponding to ~7 pc 
from the MCELS Ha image shown in Figure[2] From the column 
density A^h shell calculated in the last section, we obtain a density 
of « s heii ~ 70cirT 3 inside the shell. Since the ISM around N 158 
shows Ha emission, its temperature i s probably around 8000 K , 
typical for the warm ionised medium (McK ee & Ostrikerll 19771) . 
resulting in a speed of sound of ~ 9.2 km/s. With the inferred 
pressure of P/k = 1.1 x /~ I/2 x 10 6 cirT 3 K from equation [3] 
we obtain a density of the ambient medium of ~13 cirT 3 , yield- 
ing a compression factor of the shock of ~6. Since the cool- 
ing time behind the shock is very small due to the high density 
in the shell, the shock will be isothermal and, as shown below, 
also strong. Since we now know the radius, energy content, and 
density, we can use the solution for a wind blown bubble by 
IWeaver et alj ( 1 1977b to find out the age of the bubble, which is 
about t 1 . 1 Myr resulting in an energy input rate over this time 
interval of L sup erbubbie = 2.3 x 10 38 erg s~ l . In addition, the mass 
inside a homogeneous bubble is M — 2.3 In x/iX mn x V — 770 
M , where fi = 0.61 is the mean molecular weight of a fully 
ionised gas and Mh is the hydrogen mass. Thus, a mass loss rate 
of 6.9 x 10~ 4 M /yr over 1 . 1 Myr is derived. In the following we 
will discuss if massive stars can account for such a large mass 
loss and energy input rate. 

In order to calculate the mass loss and energy input rates 
of OB stars (in cgs units) we use t he mass-luminosit y relation 
for stars with 10 < M/M Q < 50 by IVitrichenko et al.1 (120071) to 
obtain the mass M and luminosity L+ of a star with bolometric 
magnitude Mb f. 



M bol = 1.6-6.91og(M/M ) 



L* = 19(M/M ) 2 - 76 L, 



The radius of the star is obtained from 



R 



Ik; 



(4) 
(5) 



(6) 



where T e g is the effective temperature of the star and cr is the 
Stefa n-Boltzmann constant. M \, \ and are taken from Table 
2 of iTestor & Niemelal (Tl998). The wind velocity is determined 



accor ding to the theory of radiation-driven winds (ICastor et al.l 
11975b : 



2GM 



R 



X(l -U/Ledd) 



0.5 



(7) 



where a 2.5 (Lamer s et al.lll995l) . The parameter v esc is the 
photospheric escape velocity and the Eddington luminosity is 
L e dd = 4nG xMxm p xc/cr T with m p being the mass of a proton 
and cr T the Thomson cross-section for the electron. The mass 
loss is determined from the single-scattering limit 



M : 



V'oo X C 



(8) 



We corrected for t he LMC metallicity (lmc following 
Leitherer et al ] (1 19921) and obtain M oc for the mass 

loss and Voo oc for the wind velocity (for hot stars with 

M > 15M ). Thus, the 67 O- and B-stars in LH 104 produce 
a mass loss of ~ 39 M or an energy input rate of Lob = 
3.4 x 10 37 erg s~ l over 1.1 Myr. 

Additionally, the WR-stars generate ~ 35 M Q or Lwr = 1.8x 
10 38 erg but only for 2.5 xlO 5 yr = 0.23 t, assuming that they 
already went through hal f of their WR-lifetime of ~ 5 x 10 s yr 
dMaeder & Mevneg|1987l) . Mass loss and wind v elocity of WR- 
binaries are adopted from iLeitherer et ail (Q997) and corrected 
fo r mass losses of WC and WN types in the LMC according 
to ICrowfherl d2007l) . Furthermore, we have to calculate the con- 
tribution of winds from WR-binaries before entering the WR- 
phase, i.e., for the remaining 0.85 Myr. With mass losses for 06 
and Q7-stars and Q4 -stars as WR-progenit ors (PR) taken from 
iGarmanv et al.l d 198 lb and velocities from Lam ers & Leithereil 
(119931) . we get 4 M Q or L PR = 8.0 x 10 36 erg s~ l for 0.85 Myr s 
0.77 t from these stars. 

In total, we find that winds can account for ~ 78 M or a me- 
chanical luminosity of L — 8.0 x 10 37 erg s" 1 over 1.1 Myr. Still, 
the largest part of the mass (~690 M ) or L SU perbubbie — L = 1.5 x 
10 38 erg s _1 corresponding to a thermal energy of 2.3 x 10 51 erg 
is 'missing', but this can be partly explained after applying an 
IMF to the star cluster. According t o the H ertzsprung Russell di- 
agram (HRD) of Testor & Niemela (1998, Fig. 6b) and using th e 
masses derived from the M-L relation (iVitrichenko et al.l l2007). 
we find that 16 stars in LH 104 have masses above 25 M in- 
cluding all O-stars, WR-binaries, the B0V stars with Id 4-41 and 
4-55, and Sk-69 259. A s an upper mass lim it we take 65 M as 
a rough estimate, since iMassev et alj (2000) suggested that the 
progenitor masses of WR-stars in this cluster should be in excess 
of 60 M @ . On th e other hand, by looking at the HRD of LH 104 
in Massey et al. (120001 Fig. 7) gives at least 20 stars in the 25 - 
65 M interval, for which spectral types or photometry are avail- 
able. Comparing both HR-diagrams shows that the total number 
of member stars is quite uncertain, but there must be at least 70 
stars above 8 M . 

Assuming that all stars formed from the same parental cl oud 
and using an IMF with r = -1.05 dTestor & Niemelalfl998l) we 
obtain a total of 74 stars between 8 and 65 M in the first case 
(16 stars with 25 < M/M e < 65) and 93 stars in the other case 
(20 stars with 25 < M/M Q < 65). In both cases there are about 2 
stars in the 65-80 M mass interval or about 3 stars in the 65 - 
90 M Q mass interval. This suggests that 2-3 supernovae (SNe) 
already exploded around 1 Myr ago and with £sn = 10 51 erg per 
SN explosion, they can easily account for the required amount 
of energy. These 2-3 SNe yield approximately 150 - 240 M of 
ejecta mass, thus 230- 3 10 of 770 solar masses in the bubble can 
be explained. Remaining discrepancies should be due to mass 
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loading and evaporation of entrained interstellar clouds and/or 
turbulent mixing of material from the cold shell. 

With the parameters determined so far we can further inves- 
tigate the geometry of the bubble and the ambient ISM. We use 
an analytical model for the expansion of a wind-blown bubble 
in an exponentially stratified medium symmetric to the galac- 
tic midplane (Baumg artner & Breit s chwer dt, in prep.) based on 
the approximation of iKompaneetsI dl960l) . With an energy in- 
put rate of L SU pertmbbie - 2.3 x 10 38 erg s 4 the bubble should 
reach a radius, i.e., semi-major axis, of 56 pc after ~1.1 Myr. 
We calculated models with different scale heights (H = 50, 100, 
and 500 pc) and obtain a density of the ambient medium of 
no = 9,11, and 13 cm 4 , respectively. Values for the height 
of the bubble above/below the galactic plane are c = 80,66, 
and 58 pc. Using b — 42 pc, we obtain a volume of the bub- 
ble of V = 2.4 x 10 51 cm 3 , 1.9 x 10 61 cm 3 and 1.7 x 10 61 cm 3 . 
Since a larger elongation, i.e., lower scale height, yields a larger 
volume, we argue that a larger scale height fits better, other- 
wise resulting in a density inside the bubble that is too low. 
With V = 1.7 x 10 61 cm 3 and norm from Sect. 13.11 we derive 
n = 3.9 x 10 4 cm 4 . This is close to the value that was used 
to calculate the mass inside the bubble and the energy input rate 
in the beginning, whereas a volume of V = 2.4 x 10 51 cm 3 for 
H — 50 pc yields only n — 3.2 x 10 4 cm 4 . The shell thickness 
in the case of H — 500 pc is ~ 5 pc resulting in a density of the 
shell of 100 cm 4 , which is somewhat higher than the density 
of 70 cm 4 derived from the observation. Finally, we calculated 
the velocity of the outer shock, which propagates into the ambi- 
ent warm medium at the top/bottom of the bubble of 31 km s 4 , 
which agrees very well with the expan sion velocity found with 
the wind solution of Wea ver et a and gives a Mach num- 

ber of M ~ 3.4. For comparison, Dunn e et al.ld2001|) reported an 
expansion velocity of the superbubble of v exp « 45 km s 4 as- 
suming a pressure-driven bubble in an homogeneous medium. 
The magn etic fields in the IS M of the LMC are on the order of 
~ 1/zG dGaensler et al.ll2005l) corresponding to an Alfven veloc- 
ity of ~0.5 km s 4 . Therefore, the effects of the magnetic fields 
are negligible. 

4. Summary 

We have found significant diffuse X-ray emission in the field of 
view of the XMM-Newton observations of the pulsar B 0540- 
69, which is not related to the pulsar. The analysis of the spec- 
trum of the emission has shown that it is purely thermal and can 
be best modelled with a hot shocked gas. The extended diffuse 
emission coincides spatially with the superbubble in the H n re- 
gion N 158 in the LMC. Therefore, we conclude that the origin 
of the X-ray emission is the hot gas inside the superbubble in 
N158. 

From the parameters obtained from the analysis of the 
XMM-Newton spectra we derive a temperature of kT = 1 keV, 
a density of n = 0.04 cm 4 , and a pressure of P/k = 10 6 cm 4 K 
inside the superbubble. These results enable us to perform ana- 
lytic calculations of the evolution of the superbubble. Since the 
OB association LH 104 t hat is located inside the superbubble has 
been studied in detail dTestor & Niemelall 19981) . the stellar pop- 
ulation inside the superbubble is well known. Using its IMF and 
HRD we calculate the mass loss and energy input rates of the 
stars. We estimate an age of ~ 1 Myr and a total energy input 
rate of L — 2.3 x 10 38 erg s 4 . The massive stars including WR- 
stars and binaries in LH 104 account for L — 8 x 10 37 erg s 4 . 
Therefore, in order to reproduce the observations, there should 
have been 2-3 SN explosions in the past 1 Myr. We have 



also performed calculations of the expansion of the superbub- 
ble for different scale heights. We show that a large scale height 
of H — 500 pc can well reproduce the observed density inside 
the superbubble as well as the density in the shell around it. The 
corresponding expansion velocity of the superbubble is also in 
good agreement with the model o f lWeaveretal](fl977l) . 
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